Eccentric black hole-neutron star mergers: effects of black hole spin and equation of 

state 
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There is a high level of interest in black hole-neutron star binaries, not only because their mergers 
may be detected by gravitational wave observatories in the coming years, but also because of the 
possibility that they could explain a class of short duration gamma-ray bursts. We study black hole- 
neutron star mergers that occur with high eccentricity as may arise from dynamical capture in dense 
stellar regions such as nuclear or globular clusters. We perform general relativistic simulations of 
binaries with a range of impact parameters, three different initial black hole spins (zero, aligned and 
antialigned with the orbital angular momentum), and neutron stars with three different equations 
of state. We find a rich diversity across these parameters in the resulting gravitational wave signals 
and matter dynamics, which should also be reflected in the consequent electromagnetic emission. 
Before tidal disruption, the gravitational wave emission is significantly larger than perturbative 
predictions suggest for periapsis distances close to effective innermost stable separations, exhibiting 
features reflecting the zoom-whirl dynamics of the orbit there. Guided by the simulations, we 
develop a simple model for the change in orbital parameters of the binary during close-encounters. 
Depending upon the initial parameters of the system, we flnd that mass transfer during nonmerging 
close-encounters can range from essentially zero to a sizable fraction of the initial neutron star 
mass. The same holds for the amount of material outside the black hole post-merger, and in some 
cases roughly half of this material is estimated to be unbound. We also see that nonmerging close- 
encounters generically excite large oscillations in the neutron star that are qualitatively consistent 
with /-modes. 



I. INTRODUCTION 

It is important to understand mergers of black holes 
(BHs) and neutron stars (NSs) not only because they are 
a chief source for ground-based gravitational wave (GW) 
detectors (such as LIGO [l|) but also because they may 
be accompanied by a diverse range of electromagnetic 
(EM) and neutrino emission. One of the more interest- 
ing possibilities in this regard is that BH-NS mergers 
may be progenitors of a fraction of observed short-hard 
gamma-ray bursts (sGRBs) 0-11 . Furthermore, existing 
and planned wide-field survey telescopes such as PTE Q , 
Pan-STARRS Q, and LSST are beginning to observe 
and classify fainter EM transient events, some of which 
are expected to be produced in BH-NS mergers from a 
variety of mechanisms (see Q for a recent, detailed ex- 
ploration of the possibilities). Coincident observation of 
these EM events with a gravitational wave signal from a 
binary coalescence would provide a wealth of additional 
information about the system beyond any individual ob- 
servation, even, for example, providing an independent 
way to measure cosmological parameters However 
this depends on having a good understanding of how to 
map the particular observations to the underlying astro- 
physical processes governing the merger. 

The great diversity in sGRBs hints at the possibility 
that there will be a corresponding diversity in the associ- 
ated GW and EM signals, and motivates the exploration 
of all viable channels for binary compact (BCO) mergers. 
BCOs may form through the evolution of primordial bi- 
naries or through dynamical capture in dense stellar sys- 
tems, such as nuclear or globular clusters. BCOs formed 



through the latter channel, which could merge with non- 
negligible eccentricity, are the focus of this work, which 
is a follow-up of a first study of dynamical capture BH- 
NS evolution using full general relativistic hydrodynam- 
ics Before delving into the details, we briefly sum- 
marize the motivation for studying this class of BCO; 
see jlOl and a related study on eccentric NS-NS merg- 
ers for more discussion. In contrast, most studies of 
BH-NS mergers to date have focused on the quasicircular 
inspiral case; see for a review of these efforts. 

A chance close-encounter of a BH and NS in a stellar 
cluster could result in a bound system if the energy loss 
due to GW emission is sufficiently large. (Tidal inter- 
action with the NS is also a source of energy loss [isj . 
though due to the relative scalings with distance, to 
leading-order the total cross section for capture is dom- 
inated by the GW emission.) Because of gravitational 
focusing, a sizable fraction of binding encounters will re- 
sult in highly eccentric binaries that merge within a few 
orbits of the initial encounter. For Newtonian hyperbolic 
orbits of systems with total mass M and relative velocity 
w at infinity, the pericenter separation rp is related to the 
impact parameter at infinity 6 by rp = lP"uP' /2M -\-0{w^) 
(unless otherwise stated we employ geometric units with 
G = c = 1 throughout). Thus, since the capture rate is 
linearly proportional to the cross sectional area, it is lin- 
early proportional to r^. For the 4:1 BH to NS mass ratio 
systems studied here (with w = 1000 km/s), estimates of 
the energy lost to GW based on the work of [31 shows 
binding encounters occur at rp < 40 M (and scales with 
mass ratio as /x^^^). The results presented here show that 
cases with initial rp < 7M (for the canonical nonspinning 
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BH), or close to 20% of binding encounters (of a 4:1 sys- 
tem with w = 1000 km/s), result in "direct collisions," 
i.e. merging on the first close-encounter. 

At first glance it may seem that such encounters would 
be too rare to be of any relevance as possible GW or 
EM transient sources; however, recent studies have sug- 
gested otherwise. In [31 estimates of BH-BH encoun- 
ters in galactic nuclear clusters suggest Advanced LIGO 
rates of 0(1 — lOf) per year, while in [l^ (following ear- 
lier work in |17Hl9l |) it was argued that direct capture 
NS-NS mergers in globular clusters could account for a 
sizable fraction of observed sGRBs. In this latter work, 
a procedure is also provided for scaling the NS-NS re- 
sults to BH-NS systems, which yields BH-NS dynami- 
cal capture event rates of 0(10 — 10^) yr^^ Gpc^'^ (us- 
ing an estimate of the relative fraction of BH to NSs 
/bh//ns ~ 0.28 as suggested in [ill). comparison, 
population synthesis models [g^] find primordial BH-NS 
merger rates from ^ 0.1 yr ~^Gpc~^ (pessimistic) to 
~ 120 yr ~^Gpc~^ (optimistic). To estimate Advanced 
LIGO detection rates would require full templates for 
these events, which we relegate to a future study; how- 
ever, in [U the signal-to-noise ratio was computed using 
post-Newtonian-based models of the early stages of the 
mergers, which showed that a subset (depending upon 
the component masses) could be observed out to 200-300 
Mpc for an average orientation, even excluding the fi- 
nal stages of the merger in the templates. This suggests 
detection rates of ^ 0.3 — 10/yr. For comparison to pri- 
mordial BH/NS binaries, [22| quotes an Advanced LIGO 
optimal detection distance of 927 Mpc for a 10 M0 BH 
to 1.4 M0 NS merger, which scales to ^ 400Mpc when 
averaging over orientation and sky location; the corre- 
sponding detection rates range from 0.2/yr to 300/yr for 
pessimistic to optimistic source population estimates. 

Though far from conclusive, there is also observational 
evidence for multiple sGRB progenitors, which could, in 
part, be due to dynamical capture vs primordial merg- 
ers. Of sGRBs with identified host galaxies, ^ 25% have 
offsets of > 15 kpc from their hosts [2^. This subset of 
sGRBs with large offsets would be consistent with kicked, 
primordially formed BCOs or with dynamically formed 
binaries in globular clusters. The latter may be preferred 
for the largest offsets 1241 . especially if primordial BCOs 
experience weak kicks [2^. Analysis of x-ray afterglows 
observed by Swift/X-Ray Telescope suggests that differ- 
ent progenitors may be responsible for sGRBs with and 
without extended emission [2^; again, one possible ex- 
planation is dynamical capture (with extended tidal tails 
leading to long-term emission) vs primordial. There has 
also been a claim that a very high-energy gamma-ray 
source observed in Terzan 5 may, in fact, be the remnant 
of a BCO merger-powered sGRB [g^l ; if true, this would 
support the claim that dense cluster environments can 
be significant sources of BCO mergers. 

It has also been suggested that some merging NS-NS or 
BH-NS binaries are actually part of a coeval or dynami- 
cally formed triple [28| . The tertiary (for example a white 



dwarf) drives the binary to high eccentricities through 
a Kozai resonance, resulting in a faster merger time. 
Though dynamically different from direct capture two- 
body systems, in Kozai-accelerated evolution the merger 
itself could take place with significant eccentricity [29[ 
and have comparable behavior at late times to the sys- 
tems studied here. However, we are unaware of any sys- 
tematic studies of expected populations and correspond- 
ing LIGO-event rates for such Kozai triple systems. 

Merging with moderate to high eccentricity could have 
significant effects on all the GW and EM observables 
from the event. In contrast to a low-eccentricity inspiral, 
the GWs arc emitted primarily around periapsis, result- 
ing in waveforms that resemble a sequence of bursts more 
than a continuous signal. Consequently, for mass ratios 
relevant to stellar mass BH-NS mergers, the evolution of 
effective orbital parameters describing the binary is not 
adiabatic, i.e. it does not occur quasistatically as in the 
early inspiral of low-eccentricity binaries. Regarding ef- 
fects associated with tidal disruption of the NS, an inter- 
esting coincidence for this class of BH-NS systems is that 
the typical radii inside of which tidal stripping begins, de- 
pending upon the NS compactness and, hence, equation 
of state (EOS), roughly coincides with the range of peri- 
center separations where the orbit becomes unstable due 
to general relativistic effects (rp < lOM depending upon 
on the spin of the black hole where M is the total mass of 
the system). The resultant dynamics could thus be very 
different from Newtonian expectations. Of course, these 
two zones also roughly correspond for quasicircular inspi- 
ral. However, highly eccentric binaries have significantly 
more angular momentum (at a given orbital separation) 
than quasicircular binaries. This additional angular mo- 
mentum could strongly affect the matter dynamics rel- 
ative to the quasicircular case, for example, resulting in 
comparatively massive disks and/or in multiple and pro- 
longed episodes of mass transfer and ejection. Ejected 
mass will decompress and form heavy elements through 
the r-process |30l - {32 |. and thus these systems could ac- 
count for a significant fraction of such elements in the 
Universe. Furthermore, subsequent decay of the more 
radioactive isotopes could lead to observable EM coun- 
terparts [1,111]. 

In [l3l we presented the first simulations of BH-NS hy- 
perbolic capture and merger for a nonspinning BH and 
a single equation of state, modeled using general rela- 
tivistic hydrodynamics (GRHD) (see [l^ for an earlier 
related study using a Newtonian-based hydrodynamics 
code) . We found a striking dependence of the outcome — 
disk mass, unbound material, GW signal — on the impact 
parameter. However, astrophysical BHs are expected to 
form with a range of spins, (see e.g. [13, [1E|) and in the 
quasicircular case spin was found to be crucial in ob- 
taining significant accretion disks (36l - [39| . Additionally, 
there is significant uncertainty about the NS equation of 
state. Since the EOS determines the NS compaction and, 
hence, the point at which the NS will become tidally dis- 
rupted, it is also an important determinant of the merger 
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outcome. Here, we expand upon the previous analy- 
sis and consider initial BH-spins aligned and antialigned 
with the orbital angular momentum (with dimensionless 
spin parameter a = 0.5 and a = —0.5, respectively), as 
well as different EOSs (the "2H" , "HB" and "B" models 
from 01 ) . Note that there is no a priori reason to expect 
alignment of the BH spin with orbital angular momen- 
tum in dynamical-capture binaries; this particular choice 
of spin direction and magnitude was purely motivated as 
a first, simple exploration of the effects of spin on the 
merger. Certainly in future studies a broader expanse 
of parameters will need to be considered. There is also 
much room for improvement with the matter description, 
including more realistic EOSs and additional physics be- 
yond GRHD. For quasicircular BH-NS inspiral [4l|, EOS 
effects were studied in (42| - |45| . the effects of magnetic 
fields i n |46| . , and higher mass ratio systems (up to 
7:1) in 

In the remainder of the paper we begin with a brief 
review of our numerical methods (Sec|IT|; more details 
are given in a companion paper [49[. In Sec. IIII[ we 
describe the particular cases we study. We focus on sys- 
tems with small initial periapsis r^, in part for the prac- 
tical reason that these binaries merge quickly and are 
thus computationally tractable, this is the regime where 
full general-relativistic (GR) effects will be most strongly 
manifested, and because maximum complementary infor- 
mation to post-Newtonian studies (e.g. [2l| or [lg|) can 
be obtained. We present the results of the simulations in 
Sec. lIVI We find that the rich variability in the dynamics 
and merger outcome as a function of impact parameter is 
compounded by considering different EOSs and values of 
BH spin. For example, we find that with prograde spin 
or a stiffer EOS significant episodes of mass transfer may 
occur during nonmerger close-encounters. Systems where 
the BH has retrograde spin or that are somewhat less ec- 
centric (than parabolic) can undergo sustained whirling 
phases before merger that are evident in the GW signal. 
As in the eccentric binary NS mergers studied in , we 
also find that strong /-mode oscillations in the NS can 
be excited in close, nonmerger encounters (see also an 
earlier study of a head-on BH-NS collision [501, though 
here the presence of the f-mode is largely a consequence 
of the initial data). In Sec. |V]we use the results of the 
simulations to calibrate a simple model for the nonadia- 
batic evolution of effective orbital parameters during each 
close-encounter. This will be used in a future study to 
build model GW templates and explore the detectability 
of these systems with gravitational wave observations. 
We conclude in Sec. IVIl 



II. NUMERICAL METHODS 

In this section, we briefly outline our numerical meth- 
ods for solving the GRHD equations. More details are 
presented in [49|, including results from an extensive 
range of tests. 



A. Evolution 

We model the evolution of BH-NS binaries by solving 
the Einstein field equations coupled to a perfect fluid. 
We solve the field equations in the generalized harmonic 
form jsil - lssi with constraint damping [s^, [l^, where 
the coordinate degrees of freedom are specified through 
source functions We employ fourth order accurate 
finite difference (FD) techniques and Runge-Kutta time- 
stepping. In order to avoid singularities we excise within 
any apparent horizons. In evolving eccentric binaries, we 
find that a damped harmonic gauge similar to the one 
described in [H^ (see also [l^l) is beneficial for achieving 
stable evolutions through merger. Specifically a damped 
harmonic gauge takes the form i/" = C("-° ~ n'^), where 
ri" is the four-vector normal to the constant coordinate 
time slices, n" is another timelike unit vector, and ^ is a 
constant controlling the magnitude of the damping. The 
|) articular form for that we found to work well is from 

We use a value of « 0.2/M, though include a spa- 
tial dependency so that ^ goes to zero at spatial infinity. 
In addition, we begin with initial data in the harmonic 
gauge and transition to this damped harmonic gauge be- 
fore the two objects begin to strongly interact. The use 
of this gauge seems to smooth out sharp features in the 
lapse that develop near merger when harmonic gauge is 
used, however, we have not studied the pure harmonic 
case in sufficient detail to conclude whether or not har- 
monic time slicing is developing a coordinate singularity. 

We describe the neutron star material as a perfect 
fluid. The fluid equations are written in conservation-law 
form and solved using high-resolution shock-capturing 
schemes. Though we have implemented several methods 
for calculating intercell fluxes and for reconstructing fluid 
primitive variables at cell faces, we used HLL [H^ com- 
bined with WENO-5 [s^ for the results presented here. 
The fluid is evolved with outflow boundary conditions us- 
ing second-order Runge-Kutta time-stepping. To resolve 
the various length scales, we use Berger- and Oliger- style 
adaptive mesh reflnement (AMR) [60] . AMR boundaries 
require special treatment in conservative hydrodynam- 
ics codes, since the fluxes emerging from a fine-grid re- 
gion, for example, do not exactly match the incoming 
fiux calculated on the coarse grid due to their differing 
truncation errors. To enforce conservation, we correct 
the adjacent coarse-grid cells using the fine-grid fluxes, 
according to the method of Berger and Colella [6lj ; more 
details on our particular implementation of this method 
is given in [49| . 

Truncation error estimates are used to generate the 
AMR-level structure. All initial data were evolved with 
a flducial "medium" resolution run, where the coarsest- 
level has 128'^ cells and covers the entire domain (we use 
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a compactified coordinate system, so this includes spa- 
tial infinity). We also chose a maximum truncation error 
so that initially six additional levels of refinement (seven 
total) are generated to resolve the BH and NS. Given 
the computational expense of these simulations, we lim- 
ited the total number of levels to seven during evolution, 
and did not allow regridding on the two coarsest-levels to 
prevent the algorithm from tracking the outgoing grav- 
itational waves beyond the largest extraction sphere of 
lOOAf. For several representative cases we also ran a 
"low" and "high" resolution simulation for convergence 
testing, where on each level the low (high) resolution 
run had a mesh spacing of 64/50 (64/96) of the medium 
resolution run, and we scaled the corresponding maxi- 
mum local-truncation error threshold-parameter used by 
the AMR algorithm assuming second-order convergence. 
Note that this procedure will not generate identical hi- 
erarchies between different resolution runs (except for 
the two coarsest levels in the wave zone, which we keep 
fixed), but on average the highest- resolution grid cover- 
ing a given coordinate cell will have the above refinement 
ratios between the different runs. To give some sense of 
the smallest scales resolved by the hierarchy, before tidal 
disruption, the low (medium, high) resolution run has 
two finest-level meshes centered around the BH and NS 
of roughly 80'^ (100^, 150^) cells each, resolving the di- 
ameter of the NS with approximately 40 (50, 75) cells 
and the BH horizon diameter with roughly 70 (85, 130) 
cells. Unless otherwise noted, results will be reported for 
medium resolution with error bars (where appropriate) 
computed from convergence calculations. 



B. Initial data 

We construct initial data by superimposing a boosted 
BH with a boosted nonspinning Tolman-Oppenheimer- 
VolkofF (TOV) star solution separated hy d = 50M. 
Though this superposition docs not strictly satisfy the 
constraint equations except in the limit of infinite sepa- 
ration, we have performed tests at various separations in 
order to verify that the superposition-induced constraint 
violation is comparable to truncation error at our reso- 
lution, in particular following an initially slightly larger 
transient that propagates away (and is partly damped 
due to the use of constraint damping) on roughly the 
light-crossing time d of the binary. In Fig. [TJ we show 
the level of constraint violation (Ca = Ha — Dxa) follow- 
ing the transient for various separations and resolutions 
for the Tp = lOM case. At these resolutions, we can 
still achieve convergence of the constraints with this su- 
perposed data, and increasing the initial separation to 
d = lOOM does not significantly affect the level of con- 
straint violation at the same resolution. This implies that 
the error introduced by the superposition is on the order 
of or smaller than the numerical truncation error at this 
resolution. 

The fact that we are beginning the binary at finite 



separation as a simple superposition of boosted single- 
compact object solutions and that there is a transient 
early time constraint- violation, both effectively introduce 
systematic errors in the parameters of the binary. To 
give some idea of the possible magnitude of this error, 
the effect on the apparent horizon-mass of the BH at 
early times and the amplitude of the density oscillation 
induced in the NS are both < 2%. Though it is difficult 
to exactly quantify how this will translate to modified- 
binary parameters, we expect it to be comparable to or 
smaller than the error introduced by setting initial or- 
bital parameters at a finite distance based on a Newto- 
nian approximation (as described in the next section). 
As we report later, the truncation error in quantities of 
interest that we do measure such as energy emitted in 
GWs is pcrcentwisc larger than this, implying that solv- 
ing the constraints or attempting more accurate initial 
representations for the metric and NS fiuid distribution 
may only offer marginal improvement in the overall ac- 
curacy of the results at these resolutions. However, for 
future higher-resolution studies it would be important to 
solve the constraints and improve the model of the phys- 
ical initial conditions for the system. 

Finally, we briefiy comment that since our simula- 
tions employ compactified coordinates such that the 
outer boundaries extend to spatial infinity, the global 
[Arnowitt-Deser-Misner (ADM)] M and J should be con- 
served. In practice, however, we must evaluate these 
quantities at a finite distance, making them subject to 
gauge artifacts, some propagating outward from the cen- 
tral BH-NS region from t = 0. For t < 200M , an ex- 
traction sphere of 300 Af is free of propagating artifacts, 
hence M (J) is conserved to better than 0.3 (2.0)% for 
all cases at medium resolution. 



III. CASES 

Motivated by possible BH-NS interactions in cores of 
galactic nuclei and globular clusters, we create initial 
data for hyperbolic encounters with varying impact pa- 
rameter. These encounters are hyperbolic since the two 
bodies have nonzero kinetic energy at (effectively) infi- 
nite separation. We will take their relative velocity at 
infinity to be w = 1000 km/s since this is the expected 
magnitude of the virial velocity in the core of a nuclear 
cluster 0,112. 

In practice, this initially positive total energy is small 
compared to the kinetic energy of the encounter itself, 
so that the orbits are nearly parabolic and have eccen- 
tricities 1 + O(10~^). In this study we also restrict our 
attention to BH-NS systems with a 4:l-mass ratio (re- 
ferring to the isolated ADM masses of the BH and NS). 
The reasons for choosing this mass ratio are in part be- 
cause it is within the range of astrophysically plausible 
values given current observations of NS and candidate 
BH masses (see, for example, [H, HI]), and in part be- 
cause within this range it is also a value where we expect 
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FIG. 1: L^-norm of the constraint violation, Ca = Ha — Oxa, 
in units of 1/M from the rp = lOM case in the lOOM x lOOM- 
rogion around the center of mass in the equatorial plane {i.e. 

\J I \\Ca\\'^<Px/ J d?x ). This is shown for low, medium, and 
high resolutions for the standard initial separation between 
the BH and NS of d = 50M; the relative magnitudes at a given 
time are consistent to good approximation with second-order 
convergence. Also shown is a medium resolution run with an 
initial separation of d = lOOM. The time is shifted so that 
the point of closest approach occurs at approximately t = 
for all cases. 



to see strong tidal-disruption effects. Certainly it would 
be of interest to explore a broader range of mass ratios, 
however, due to limited computational resources we leave 
that to a future study. 

We choose the initial positions and velocities for the 
BH and the NS according to the Newtonian equations for 
a hyperbolic orbit (see, e.g., [65|) with w = 1000 km/s 
and an initial separation of d = 50M. We vary the ini- 
tial spin of the BH, considering the values a = —0.5, 0, 
and 0.5 where negative (positive) values indicate retro- 
grade (prograde) spin in relation to the orbital angular 
momentum. We also consider three different broken F- 
law model EOSs labelled 2H, HB, and B in [43|. For the 
prototypical 1.35 M© NS that we use, these EOSs give 
compactions M^s / Rns = 0.13,0.17, and 0.18, respec- 
tively. Note that while the compaction of the B EOS NS 
is only slightly smaller than that of the HB, it has a max- 
imum mass of 2.0 Mq (the 2H and HB have maximum 
masses of 2.83 and 2.12 M©, respectively). Given recent 
observations [6^, B is therefore on the soft end of the 
allowed range for the EOS family considered here. We 
include a thermal component in the EOS, a F-law with 
Fth — 1.5 to allow for shock heating. 

We do not consider large impact parameters corre- 
sponding to initial Vp > 15M. Even so, for impact pa- 
rameters at the upper end of the range we do evolve 
the eccentricity is sufficiently large (though < 1) after 



the first close-encounter that it would be very expen- 
sive to evolve to the second encounter. To help cali- 
brate our model for orbital parameter evolution that will 
be introduced in Sec. |Vl we also consider a set of runs 
with initial orbital parameters for a bound orbit with ec- 
centricity e ~ 0.75 instead of the hyperbolic orbit with 
e w 1. These simulations can be seen as correspond- 
ing to systems that have already undergone one or more 
close-encounters and evolved to these orbital parameters. 

To keep the parameter space at a manageable size, we 
vary only one of the three parameters of BH spin, NS- 
EOS, and initial eccentricity at a time from our base 
case — an initially nonspinning BH (a = 0), a NS with 
the HB EOS, and the two objects initially with New- 
tonian orbital parameters corresponding to a marginally 
unbound orbit (e « 1) — and then consider a range of 
impact parameters. See Fig. [2] for plots of the NS tra- 
jectory for several cases, and Fig. |3] for snapshots of the 
rest-mass density at select times illustrating aspects of 
the matter dynamics. 



IV. RESULTS 

Varying the parameters as discussed in the previous 
section, there is much degeneracy in the qualitative fea- 
tures that arise (which will need to be addressed in future 
studies investigating extraction of source properties from 
GW and EM observations). This is essentially because 
the leading-order source of the variability is rooted in 
the following two properties of the system: (1) the NS 
radius, varied by altering the EOS (as the NS mass is 
fixed in this study), (2) the location of the innermost 
stable orbit (ISO) varied by changing the spin of the 
BH or the eccentricity of the encounter. For equatorial 
geodesies on a black hole background, the ISOs asymp- 
tote to circular orbits, though these should not be con- 
fused with the innermost stable circular orbit (at r = 6M 
in Schwarzschild) , or ISCO. The circular orbits (in the 
range r — 3M to 6M in Schwarzschild) associated with 
the ISOs with nonzero eccentricity are unstable, and un- 
der infinitesimal perturbation their noncircular nature is 
manifest in the form of zoom- whirl behavior [67| . Specif- 
ically, depending upon the size of the perturbation, the 
geodesic will undergo a number of near circular orbits 
(the "whirls"), followed either (depending upon the sign 
of the perturbation) by a plunge into the BH or by a sin- 
gle near elliptical orbit (the "zoom"), and in the latter 
case the motion repeats. Away from the geodesic limit 
there is still an effective ISO, where qualitatively simi- 
lar behavior occurs [g^, though radiation reaction will 
eventually drive the system to a merger for all bound 
systems. 

In terms of the gravitational dynamics, the closer the 
periapsis is to the ISO, the more whirling that occurs, 
resulting in enhanced GW emission and more rapid evo- 
lution of effective orbital parameters. Regarding the mat- 
ter dynamics, the ISO is essentially the "event horizon" 
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for fluid elements following geodesies, or being close to 
geodesic. What this implies is on a close-encounter, if 
the NS is not disrupted and the periapsis of the orbit is 
within the ISO a merger will result and the entire NS will 
fall into the BH. If the NS is disrupted, following the es- 
sentially Newtonian redistribution of angular momentum 
that results, fluid elements within an ISO corresponding 
to their effective eccentricity will immediately plunge into 
the BH, while the rest will either move out onto eccen- 
tric orbits to later fall back onto an accretion disk or 
be ejected from the system. If the majority of the NS 
mass ends up outside the ISO following disruption, once 
it moves beyond the tidal-disruption zone there will be 
sufficient self-gravity for the material to recoalesce into 
a NS core. Similarly, in a partial disruption where only 
the outer layers are stripped from the NS, some stripped 
material will accrete into the BH, some back onto the 
NS, and a portion will be ffung out unbound. 

In the following subsections, we break down the dis- 
cussion of phenomenology of the encounters by first 
summarizing the results from the base case presented 
in (Sec. IIVX)) . then describe in turn what happens 
when sp in (Se c. IIVB|) . EOS (Sec. IIVC|) . and eccentric- 
ity (Sec. IIVD[) are changed relative to the base case. In 
order to have some intuitive understanding of why the 
changes have the effects they do, it is useful to keep the 
above discussion in mind. This can also allow one to 
anticipate what would happen if parameters we are not 
altering here are changed, for example, the BH or NS 
mass. 



the burst of GWs associated with the first periapsis pas- 
sage. For one case in each class {vp = 5.0, 7.5, 10.0) we 
performed a convergence study which showed approxi- 
mately second-order convergence and allowed us to per- 
form a Richard extrapolation to estimate errors in the 
resulting GWs. To give further indication of the trun- 
cation error in these runs, in Fig. 2] we plot the error in 
the trajectory of the NS in a fiy-by case = 10 (see 
Fig. [2] for the medium resolution trajectory). The GW 
energy and angular momentum emitted (including ex- 
trapolated values from the resolution studies) as well as 
the disk properties of those cases followed through merger 
are summarized in Table ID Table |ll] shows the spin of 
the post-merger BH (for this base set of runs as well as 
the subsequent parameter survey) for the cases we fol- 
lowed through merger. As the threshold in Vp dividing 
(1) and (2) is approached there is a dramatic enhance- 
ment in the gravitational energy and angular momentum 
released during the close-encounter. 

The amount of material remaining after merger, which 
could potentially form an sGRB-powering accretion disk 
also depends significantly on impact parameter. Below 
the threshold dividing (1) and (2), there is a sizable 
amount of remaining material in excess of 20% in the 
Tp = 6.81 case. In most cases « 50% of the material is 
unbound. A simple numerical estimate of the fallback 
time for the bound material based on when the elliptic 
orbit will return to the accretion disk shows the charac- 
teristic t~^/^ scaling [69|. 



A. Zero-spin survey with HB EOS 

Here we summarize the results obtained for simula- 
tions with an initially nonspinning BH and a NS with 
the HB EOS (see [10| for more details). We considered 
a range of periapsis separations from rp/M = 5.0 to 15 
{i.e., 50 to 150 km). Henceforth, we will consider rp 
to be normalized by M. In all of these cases, sufficient 
energy is carried away by GWs to result in a bound sys- 
tem. Our simulations exhibit three types of behavior: 
(1) a direct plunge {rp = 5.0,5.83,6.67,6.81); (2) follow- 
ing the initial periapsis passage, a single elliptical orbit 
and then a plunge (r^ = 6.95, 7.22, 7.5); and (3) following 
the initial periapsis passage, a long-period elliptical orbit 
{rp = 8.75, 10.0, 12.5, 15.0). For the latter group (and the 
high-resolution rp = 7.5 run), we do not simulate the en- 
tire orbit since the length of such simulations would make 
them very computationally expensive ^ , and we focus on 



^ For example, based on an estimate using the emitted GWs and 
assuming a Newtonian orbit the Vp = 8.75 case will undergo 
another close-encounter after ~ 7000 M, which would take 2- 
3 months of wall clock time to simulate at medium resolution 
(while the larger Vp cases will take even longer). Though such 
long runs are not unheard of, we chose to use our limited com- 



B. Effects of black hole spin 

As expected from Kerr geodesies, prograde BH spin re- 
sults in a smaller critical-impact parameter below which 
the NS merges with the BH on the first encounter than in 
the nonspinning case; the converse is expected for cases 
with retrograde initial BH-spin. For the cases with initial 
BH-spin a ~ 0.5 (top of Table llll| . only the rp = 5.0 case 
merges on the first encounter while rp = 5.5, 6.0, and 6.25 
go back out on a short elliptic orbit before merging. The 
rp ~ 7.5 and 10.0 cases go out on a long elliptic orbit 
after the initial periapsis passage, and we did not follow 
these through merger. For the rp = 5 case, we also per- 
formed convergence tests in order to estimate truncation 
error for a spinning BH where the effective ISO is closer 
to the BH. From this we compute Richardson extrapo- 
lated values for the energy and angular momentum of 
the resulting GW as well as for the amount of matter left 
over post-merger; these values are reported in the table. 
In Fig. [5] we also show the amount of rest- mass exterior 
to the BH horizon as a function of time and resolution 
for this case. 



putational resources to explore a greater number of parameters. 
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FIG. 2: Trajectories of the NS center-of-mass for various simulations. The r-p = 7.0 with 2H EOS, r-p = 7.0 with B EOS, and 
r-p = 5.5 a = 0.5 (HB EOS) both undergo a close-encounter followed by a short elliptic orbit before merging. Note that while 
in both the Tp — 7.0-cases the NS approaches the BH on essentially the same orbit, the dynamics around the close-encounter 
and subsequent orbits are very different due to EOS effects (see Sec. IIV C)) . The Vp = 10 (HB EOS) undergoes a long-period 
elliptical orbit following the initial periapsis passage. The remainder of the cases shown merge on the first encounter while 
displaying various degrees of whirling behavior. 



One important consequence of the reduced critical- 
impact parameter in the prograde case is that around 
this threshold, the tidal forces on the NS arc greater and 
the resulting accretion disks consequently larger. Even 
more striking, the enhanced tidal interaction can lead to 
significant mass transfer from the NS to the BH even 
in nonmerger close-encounters, as, for example, with the 
Vp = 5.5, a = 0.5 case. Here the NS becomes highly 
distorted and loses approximately 16% of its mass to the 
BH (see Fig. |6|). Hovirever, a compact (albeit highly dis- 
torted) star remains (see Fig. |3] middle, bottom) until it 
merges with the BH on its second close-encounter. 

For the a = —0.5 cases (see bottom of Table [Hi]) . = 
5.0, 7.5, and 8.13 are direct mergers while the Vp = 8.28 
case goes through a single elliptic orbit before merging. 
The Tp = 8.44, 8.75, and 10.0 cases go out on longer ellip- 
tic orbits which we did not follow to completion. Because 
of the larger critical-impact parameter for merger on the 
first encounter, there is less tidal disruption and, hence. 



the accretion disks contain < 1% of the total NS rest- 
mass for all the cases followed through merger. However, 
this also means that the BH and NS can undergo more 
whirling behavior in the critical regime before merger. 
This is especially evident in the waveform for the nearer 
threshold merger case of Vp = 8.13, a = —0.5 as shown in 
Fig. [TUl The gravitational wave signal shows several cy- 
cles of almost constant amplitude and frequency, indica- 
tive of a nearly circular orbit. Compared to the near- 
threshold nonspinning case (r^ = 6.81, also shown in 
Fig. I10|) . this whirling period is much more pronounced 
in the waveform. 



C. Effects of equation of state 

We also consider the B EOS, which is softer than the 
HB, and the 2H EOS, which is stiffer. Though the dif- 
ferent EOSs have only a slight effect on the critical im- 
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FIG. 3; Rest-mass density in the equatorial plane from various BH-NS simulations, left to right, top to bottom: (1) the BH 
and NS undergoing a close-encounter {t = 242 M, Vp — 7.0, B EOS), (2) the NS undergoing a whirling-phase before merging 
{t = 269 M, Tp — 8.13, a — —0.5), (3) the NS stretched into a long tidal stream during merger {t — 506 M, Vp = 7.5, e = 0.75), 
(4) a mass transfer episode {t = 292 M, rp — 7.0, 2H EOS), (5) towards the end of a mass-transfer episode during the NS's first 
periapsis passage (t = 272 M, rp = 5.5, a = 0.5), and (6) a nascent accretion disk (t = 388 M, rp = 5, a — 0.5). Recall that rp 
is reported in units of total mass M. The color scale is logarithmic from 10"*^ to 1 times the initial maximum density (pmax = 
8.3 (3.7,9.8) xlO^* g cm-^* for the HB (2H, B) EOS). The BH is roughly the same coordinate size (with diameter 3 M) in 
all panels, which can be used to infer the relative scale of each snapshot. 
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TABLE I: Disk properties and GW energy and angular-momentum losses for an initially hyperbolic (e ~ 1) encounter of a 
zero-spin BH and NS with HB EOS. Dashed entries correspond to cases that merge during the first encounter, and hence have 
no "first periapsis"; dotted entries correspond to binaries that were only evolved through first periapsis passage. 
" Rest-mass remaining outside the BH shortly 50 M) after the end of merger, normalized by the initial total rest-mass. 

Unbound rest-mass estimated using local fluid velocities and assuming a stationary metric. 

Rough initial accretion timescale (race = Mo/ Mq) evaluated shortly after merger. 

Energy and angular momentum lost to GWs during the first close-encounter. 

Total GW energy and angular-momentum losses for cases which were followed through merger. 

Results are from medium-resolution runs; values in parentheses are Richardson-extrapolated estimates using low and high 
resolutions, where available. Note that the relatively large error for rp = 7.5 (and to a lesser extent rp = 5, 10) is due in part 
to truncation error altering the actual periapsis by a small amount, and in this regime the GW emission is highly sensitive to 
binary separation (Fig. I12p . 

^ For the = 5 case Mo/Mo{t = 0) was the same at the three resolutions to within ~ 5 x 10"^; however, for such low disk 
masses we expect systematic effects, e.g. the numerical atmosphere to be important. 
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FIG. 4: Convergence of the NS trajectory for the rp — 
10, a = 0-case (see Sec. IIV A[) . What is plotted is the mag- 
nitude of the difference in the center of rest-mass position of 
the NS between the low and mediumi-resolution runs, and be- 
tween the medium and high-resolution runs scaled assuming 
second-order convergence. The close-encounter occurs near 
t — 230Af , where the largest "perturbation" of the orbit due 
to truncation error happens; the relatively large growth of the 
difference in trajectories following this is mostly just a reflec- 
tion of this change in the orbital parameters that occurred 
around periapsis. 



a = HB EOS 



rp 5.00 6.67 6.81 6.95 7.22 7.50 
aflnai 0.49" 0.45 0.37 0.47 0.50 0.50 
a^s 0.40 0.46 0.47 0.47 0.48 0.49 



a = 0.5 


a = -0.5 


rp 5.00 5.50 6.00 6.25 
afinai 0.74 0.71 0.71 0.70 
aefi 0.72 0.74 0.77 0.77 


5.00 7.50 8.13 8.28 
0.16 0.25 0.22 0.24 
0.08 0.17 0.19 0.19 


2H EOS 


B EOS 


rp 5.00 6.75 7.00 
ciflnai 0.50 0.29 0.33 
flefi 0.40 0.46 0.47 


5.00 6.25 7.00 
0.48 0.52 0.48 
0.40 0.45 0.47 


e = 0.75 
rp 7.50 7.81 
afinai 0.44 0.49 
aofi 0.46 0.47 



"afinai = 0.49 ± 0.01 

TABLE 11: Post-merger BH spin (dimensionless) for various 
initial conditions. Also shown is the effective spin a^fi, used 
in the model described in Sec. |V] and calculated from initial 
conditions using (j3}. 



FIG. 5: Effect of resolution on the disruption and subsequent 
accretion of NS material, as measured by the total rest-mass 
exterior to the BH horizon, for the = 5, a = 0.5 case (see 
Sec. IIV B|) . The amount of rest-mass remaining at late times 
at the different resolutions is consistent with approximately 
second-order convergence. 
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FIG. 6: Top: The amount of rest-mass (normalized to the 
total rest-mass of the NS) accreted by the BH as a function 
of time for the rp — 7.0 2H-case and for the rp = 5.5, a — 0.5 
HB-case. In both cases there are two significant episodes of 
mass transfer. Bottom: The angular momentum of the BH 
horizon in units of the total ADM-mass squared, , as a 
function of time for the same two cases. 
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TABLE III: Disk properties and GW energy and angular-momentum losses for an initially hyperbolic (e « 1) encounter of a 
BH with spin a = ±0.5 and NS with HB EOS. The same comments and set of footnotes " to ^ apply as in Table |T] 
^ For the rp = 5, a — +0.5 case the Richardson extrapolated value of Mo,u/Mo{t = 0) was computed using just the medium 
and high-resolution results and assuming second-order convergence due to the low amount of unbound material in the 
low-resolution case. 
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FIG. 7: The fraction of rest-mass outside a given radius 
(1/Afo) pdV for an isolated star with the HB or 2H EOS 
and M = 1.35 Mq. Here radius r is normalized to the radius 
of the HB star (-Rhb). The horizontal, dotted-lines indicate 
the amount of material accreted into the BH during the initial 
close-encounter for the rp = 5.5-case with HB EOS (bottom) 
and the Vp — 7.0-case with 2H EOS (top) as shown in Fig. [S] 



pact parameter, they can have a dramatic effect on the 
near-threshold behavior. For the B EOS, we chose four 
different impact parameters (bottom of Table lIVp : rp = 
5.0 and 6.25 (direct mergers); rp = 7.0 (after periapsis 
passage, a short elliptic orbit and then merger); and rp = 
7.5 which we did not follow through its full elliptic orbit 
after the first close-encounter. The more compact B EOS 



NS experiences less tidal disruption and forms accretion 
disks with < 1% of the total NS rest- mass for all cases fol- 
lowed through merger. For the 2H EOS, we considered 
five different impact parameters (see top of Table IIVP : 
direct mergers rp = 5.0 and 6.75; rp = 7.0 (after periap- 
sis passage, a short elliptic orbit and then merger); and 
rp = 7.25 and 7.5 which we did not follow through the 
full elliptic orbit after the first close-encounter. 

The less compact 2H EOS NS experiences significant 
tidal deformation, and for the rp = 6.75 and 7.0 cases 
there is « 30% of the NS material leftover post-merger. 
In the case of rp = 7.0 the deformation and tidally- 
induccd oscillations of the NS are especially severe (sec 
Fig. [21 bottom left) and there is a period of significant 
mass transfer (Fig. [6]) from the NS to the BH during 
the first close-encounter. This is very similar to the 
rp = 5.5, a ~ 0.5 case, though more pronounced. For 
both the 2H EOS and the a = 0.5 with HB EOS sim- 
ulations no significant mass transfer is found for runs 
with rp sHghtly above 7.0 and 5.5, respectively, so we can 
use these as an effective measure of the separation da at 
which the star begins to overflow its Roche lobe. Though 
the individual values are somewhat higher than expected 
from Newtonian theory, it predicts a scaling with com- 
paction dfi (X C^5, where Cms = Mns/Rns', this does 
approximately hold here since Chb/C2h = 1-3. This 
ignores any relativistic effects, both gauge and physical, 
the latter including, for example, the effect of BH spin 
on tidal disruption. In Fig. [7]wc show the mass profiles 
of an isolated star with the 2H and HB EOS. This shows 
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FIG. 8: NSs exhibiting large-amplitude /-mode oscillations following the first interaction with the (nonspinning) BH. The 
cases shown are rp = 6.95 with the HB EOS at t = 341M (left) and r-p = 7.0 with the 2H EOS at t = 356A4" (for the latter case, 
see the bottom left-panel of Fig. [3] for a snapshot of the density near the time of closest approach to the BH when this large 
oscillation is excited). The color map shows the rest-mass density in the equatorial plane on a logarithmic scale, normalized to 
the instantaneous maximum density. The arrows show the velocity in the NS center-of-mass frame. The longest velocity arrows 
correspond to velocity magnitudes 0.084 (left) and 0.17 (right). Though the two cases have similar periapsis radii, the lower 
compaction in the 2H-case leads to a much stronger tidal interaction and an /-mode, which shows nonlinear characteristics. 



that the outer spherical shell that contains the approxi- 
mately 40% of the NS material that is accreted into the 
BH during the first close-encounter for 2H corresponds to 
a volume containing almost no material in the HB case. 
Again, though we are ignoring the complicated details of 
mass-transfer dynamics, this is suggestive as to why, for 
the BH-NS system studied here, there is significant mass 
transfer for the 2H EOS around rp = 7.0 but very little 
for the HB EOS. 

Since for the 2H rp = 7.0 case so much mass is trans- 
ferred before the merger and there is a strong disruption 
during merger, the gravitational wave signal resulting 
from merger itself is significantly weaker than for other 
cases (sec Fig. [T0|) . The full waveform from Fig. [11] also 
shows that the GW pulse from the fiy-by dominates the 
signal compared to the merger part of the waveform. 

Between the initial close-encounter and merger there is 
also evidence in the GW signal of excited /-modes within 
the NS; this was also observed in eccentric-binary NS en- 
counters [ll| . Here we give a qualitative account of the 
stellar dynamics to show that the dominant oscillation is 
akin to that of an /-mode of an isolated, perturbed star — 
a detailed study, aside from the difficulty of applying a 
perturbative analysis in such a transient and in some 
cases highly distorted star, is beyond the scope of this 
study. Figure [5] demonstrates the density and velocity 
distortions for the rp — 7.0, 2H case (right panel) as well 



as the rp = 6.95 HB case in which the distortion is less ex- 
treme (left panel). The latter bears a particularly strong 
resemblance to the pure I = m = 2 /-mode fiow pattern 
(see, e.g., Fig. 19 of [z3|)- In both cases, an animation 
of the density field seems to suggest a rotating, distorted 
NS. However, as Fig. [5] demonstrates, the velocity pat- 
tern is not one of overall rotation but rather that of an 
oscillatory mode. Indeed, the circulation theorem should 
hold to a good approximation in the bulk of the NS mate- 
rial (though there are entropy-generating shocks near the 
surface). Thus, the tidal interaction with the BH would 
not induce rotation in the usual sense, but rather oscilla- 
tory modes with rotating patterns. Many such modes are 
likely to be excited, but the I = m = 2 /-mode seems to 
dominate. We performed a spherical-harmonic decompo- 
sition of the star's rest-mass density p = J2 C^imYim at se- 
lect times and found C22 to be the largest coefficient next 
to Coo- We also checked that the perturbation amplitude 
grows monotonically with radius. For the rp ~ 6.95 sim- 
ulation with the HB EOS IC22/C00I « 0.01 - 0.02 at 
r = 0.4A/ grows to IC22/C00I « 0.15 - 0.20 at r = 0.8M. 
Here r is the radius of the sphere centered on the NS 
center-of-mass on which the coefficients are calculated. 
The next largest coefficient C20 is smaller by a factor of 
> 2. The rp = 7.0 simulation with the 2H EOS shows 
similar behavior although the coefficients are somewhat 
larger with IC22/C00I ~ 0.03 and 0.25 at r = 0.4M and 
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0.8M, respectively. 

We briefly comment on the possible detectability of 
such an /-mode excitation in GWs. Such an observation 
in principle could provide a wealth of information about 
the structure of the NS, in particular, since /-mode fre- 
quencies are quite sensitive to the EOS (see for exam- 
ple [7l|). However, here, (a) the amplitudes are quite 
low relative to the dominant GW emission (see Fig. [TT|) , 
(b) for the cases where the largest amplitudes are excited, 
the initial rp is sufficiently small that only a few cycles of 
waves will be emitted before subsequent merger, limiting 
the signal-to-noise that could be built up, and (c) the 
frequency is quite high (above 1 kHz) and thus not in a 
regime where the AdLIGO-class detectors are very sensi- 
tive. Thus, even if there is a sizable population of eccen- 
tric merger events as studied here, it is unlikely that any 
corresponding /-mode excitation will be observed with 
the current generation of GW detectors. 



D. Bound eccentric evolution 

Finally, we simulated binaries where the BH and NS 
have initial-orbital parameters corresponding to a bound 
orbit with e = 0.75 for a range of values of r^. These 
simulations can be viewed as corresponding to systems 
that have already become bound through one or more 
close-encounters and lost some of their initial eccentric- 
ity. In practice, the length of such evolutions would be 
very computationally expensive to follow in full. Of the 
impact parameters considered (see Table IVT) . rp = 7.5 is 
the only direct-plunge, and rp = 7.81 is the only case 
we followed past periapsis passage to merge on the sec- 
ond close-encounter. For the remaining cases {rp = 8.13, 
8.75, and 10.0), we only followed partially through their 
elliptic orbits after the first close-encounter. The thresh- 
old for merger on the first encounter moves out slightly in 
this case as an appeal to geodesies in Schwarzschild would 
suggest. In Fig. [HI we plot the gravitational waveforms 
from two fly-by close-encounters Vp = 8.13 and 10.0. In 
addition, we fit the expected waveform according to a 
Newtonian-order quadrupole approximation [3, IH, [ill 
(NQA) to our numerical results by multiplying by an 
overall factor. As found in [loj and in our other simula- 
tions here for the marginally unbound case, though the 
shape of the waveforms agree quite well, the numerical re- 
sults exhibit a significant amplitude enhancement that is 
larger the closer one gets to the threshold for merger dur- 
ing the close-encounter. In addition, from Fig.[TU]we can 
see that the near-threshold merger waveform shows more 
evidence of whirling than the initially unbound case. 



V. EVOLUTION OF ORBITAL PARAMETERS 

As noted in [loj . there is a dramatic enhancement in 
gravitational-wave energy and angular momentum emit- 
ted in a BH-NS close-encounter that occurs as the thresh- 
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FIG. 9: The real and imaginary components (black diamonds 
and red squares) of the / = 2, m = 2 spherical-harmonic of 
r*4 for e = 0.75 and Tp = 8.125 (top) and Tp = 10 (bottom). 
For comparison the NQA analytical results are shown multi- 
plied by an overall factor so that the magnitude and phase 
match at peak {t — 0). 



old for merger on the first encounter is approached. In 
this section, wc will attempt to explain this enhancement 
by analogy to the zoom-whirl orbits in Kerr spacetimes, 
and extrapolate this behavior to close-encounters with 
different eccentricities and spins. Then we will use these 
results to approximate the evolution of orbital param- 
eters for BH-NS systems undergoing a scries of close- 
encounters. 



A. Zoom-whirl enhancement 

As discussed in the previous section and as can be seen 
in Fig. [121 though different values of initial spin or eccen- 
tricity may shift the critical value of Vp, which we shall 
call Tc at which the NS merges with the BH during a 
first encounter, for all the cases considered there is a sig- 
nificant enhancement in the energy of the gravitational 
waves emitted as this threshold is approached. This can 
be understood by analogy to the zoom-whirl orbits in 
Kerr spacetimes. As the BH-NS approaches the critical- 
impact parameter, it is closer and closer to the unsta- 
ble orbit dividing plunging and nonplunging orbits and 
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TABLE IV: Disk properties and GW energy and angular-momentum losses for an initially hyperbolic (e « 1) encounter of a 
BH with zero-spin and a NS with 2H and B EOS. The same comments and set of footnotes " to ■'^ apply as in Table U 
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TABLE V: Disk properties and GW energy and angular-momentum losses for an initially eccentric e = 0.75 encounter of a 
zero-spin BH with a NS with HB EOS. The same comments and set of footnotes " to apply as in Tabled 




time (M) 

FIG. 10: The real and imaginary components (solid, black- 
lines and red, dotted-lines) of ^'4 on the z-axis (perpendicular 
to the orbital plane) during merger for the following cases (top 
to bottom): rp = 5, HB; rp = 6.81, HB; rp = 7.5, e = 0.75, 
HB; rp = 8.13, a = -0.5, HB; and rp = 7.0, 2H (see Fig. [ll] 
for the full signal). The waveforms are aligned so that the 
peak occurs at t = with zero phase. 



will, therefore, exhibit more and more whirling behavior. 
Though the effect is most dramatic near threshold, this 
zoom-whirl analog is a way to understand the significant 
enhancements over the NQA predictions that persist till 
Tp « 10 (see Fig. 4 in [l3|)- This motivates a fit to the 
GW energy emitted in a whirling close-encounter with 
the functional form (see Fig. [T^ 

Egw = Eoil - {Srp/Ar) (2) 

where Srp — rp~rc, Eq is the difference in energy between 
a quasicircular orbit and an eccentricity e orbit both with 
''p = ^c, A is the range over which zoom- whirl-like be- 
havior dominates the GW- emission energetics, and 7 is a 
parameter that in the geodesic analogue is related to the 
instability exponent of the corresponding unstable circu- 
lar orbit. Given the limited number of points we have, 
we choose to use 7 as our single fitting parameter. We 
set Tc to the average of the closest sub and supercriti- 
cal values of rp from the simulations. We choose A = 3 
by inspection so that values of rp outside this range are 
well-approximated by the NQA model. Certainly, in fu- 
ture work once more data points are available it would 
be preferable to more systematically fit to the other pa- 
rameters as well, and refine the model. 

The top panel of Fig. [1^ shows this fitting performed 
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FIG. 11: The real and imaginary components (solid black 
lines and red dotted lines) of ^'4 on the z-axis (perpendicular 
to the orbital plane) at r = 70 M for the rp = 7.0 2H simu- 
lation. The large burst at t « 300 M comes from the initial 
fly- by, where the NS becomes extremely distorted (Fig. [3]) and 
loses a significant portion of its mass (Fig. [6]). The smaller 
pulse at t = 600 M comes from the merger. In between, there 
is a smaller component of the signal coming from the tidally 
induced oscillation of the NS. 



individually to each set (e, a) of simulations run. We 
should note that due to limited computational resources 
we have not been able to perform this fitting at multiple 
resolutions in order to be able to estimate the effect of 
truncation error on 7. However, due to the sensitivity 
of Eqw on rp that can be seen in Fig. [121 one would 
expect that the truncation error in these values will be 
dominated by the resolution dependence of Tc as opposed 
to the error in 7. In the remainder of this section we will 
use the data from the e « 1, a = set to calibrate 
geodesic-inspired extensions to arbitrary values of (e, a), 
which can then be compared to the other three sets of 
simulation data (bottom panel of Fig. [T2|) to test how 
well this extrapolation works. 

To extend this model to different values of eccentricity 
e and BH spin a, we assume that the dependency of Tc 
and 7 on these values is the same as in the geodesic case 
on a Kerr background. We will use QeS to refer to the 
effective spin parameter we use in the Kerr formulas. For 
better correspondence with the Kerr spacetime, we want 
to take Oeff to be the approximate spin of the BH that 
would form if a merger occurred. We do this rather then 
use the initial spin of the BH based on results from , 
which suggest that far from the geodesic limit the total 
angular momentum of the binary is more important than 
the BH spin. We estimate Ocft using 



Ocff = ao 



-a(il/BHA/)' 



(3) 



where Mbh, o-, fp, and e are the initial BH-mass, initial 



BH-spin, periapsis, and eccentricity, respectively, of the 
encounter for which we want to compute an a^s, and 
flo ~ 0.5 is the final spin measured from the e sa 1, a = 0, 
= ''cO ~ 6.9 simulation. In Table [Hi we show how 
Coff compares to the final BH-spins in the simulations 
that we followed through merger. This simple formula 
based on the total angular momentum of the system docs 
not attempt to capture any of the complications due to 
differences in matter dynamics or gravitational radiation 
between the different cases, and as can be seen does not 
always capture the trends with Vp seen in the simulations. 
We use Oeff because it has a simple motivation, and it 
does a decent job of estimating the final BH-spin for the 
purposes of this model. This allows us to extend this 
model to other values eccentricity e and BH-spin a that 
were not simulated. 

Recall that for equatorial geodesies with eccentricity e 
in Boyer-Lindquist coordinates with BH-spin parameter 
a, there is a value for the periapsis that corresponds to 
a marginally unstable orbit r^^(e, a). It can be found 
by solving the equation (see, for example, [t^) {r^^Y = 
(J — aEY[\ + e)/(3 — e), where E and J are the orbit's 
specific energy and angular momentum, respectively. For 
our model, we assume that rc(e, Oes) oc r'^'^(e, Oeff), and 
we fit the proportionality constant using our numerical 
results for e w 1, a = 0. 

The instability exponent for unstable circular orbits in 
Boyer-Lindquist coordinates is given by [H, [t^ 



7"^(e,a) 



(4) 



r 

2^ 



4M 



ir^D H -{rR^u^ - AMauj - r + 2m) 



where i? = + 0^(1 + 2M/r), D 



-1/2 



2Afr, 



M/{Ma ± vMr^) (± for prograde and retrograde, 



respectively), and we set r 



e, a). Again, we assume 



that 7(e,aeff) oc 7^^(e,aeff) and fit the proportionality 
constant using our numerical results for e « 1, a = 0. 
For the angular momentum lost to gravitational waves 
in a close-encounter we assume a similar expression 



Jgw = Jo(l - {Srp/Ay) 



(5) 



where again Jq is the difference in angular momentum of 
a quasicircular orbit and an eccentricity e orbit evaluated 
at the same separation rp. 

This simple prescription for estimating gravitational- 
wave energy and angular-momentum loss has some ob- 
vious limitations. We are extrapolating the critical im- 
pact parameter and instability exponent based on a Kerr 
spacetime to a BH-NS spacetime that is dynamic and 
nonperturbative. We are also ignoring tidal effects or de- 
pendence on EOS in this model. Nevertheless, in Fig.lT^ 
we compare how closely this predicted scaling with spin 
and eccentricity matches that from simulation results. 
Given its simplicity, as well as the numerical error in 
these results, this model does a satisfactory job of cap- 
turing the trends in these scalings. 
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FIG. 12: Top: Energy lost to GWs during the initial close- 
encounter (i.e. excluding merger) as a function of rp for initial 
BH spin a = 0,-0.5, and 0.5 and for initial eccentricity e = 
1 and 0.75. The functional form £'o(l - {Srp/Ay) (lines) 
motivated by zoom-whirl dynamics is a fit to the simulation 
results (points). Svp = rp — Vc where rc is the threshold 
value for merger during the encounter. Eo is the difference 
in energy between a quasicircular orbit and an e ~ 1 (0.75) 
orbit, both having Vp = r^. A is the range over which zoom- 
whirl-like behavior dominates the GW-emission energetics. 7 
is a parameter that in the geodesic analogue is related to the 
instability exponent of the corresponding unstable circular 
orbit, here, we use it as our fitting parameter. We obtain 
7 = 0.19, 0.13, 0.25, and 0.16 for (e, a) = (1,0), (1,-0.5), 
(1,0.5), and (0.75,0), respectively. Bottom: This shows the 
same simulation data points as the top figure, though here we 
use the e = 1, a = case (solid, black points) to determine 
the free parameters for the method described in the text to 
extrapolate the values of and 7 to the other three cases. 



B. Systems undergoing multiple close-encounters 

The enhanced GW energy and angular-momentum 
losses during a close-encounter for a given Vp result in 
more rapid loss of eccentricity and larger rate of decrease 
of Vp of the next encounter. Figure [13] shows approx- 
imate trajectories in e and rp for binaries on initially 
marginally unbound (i.e. e = 1 + e) orbits for a range 
of initial BH spin. These results were obtained using the 
above model for energy and angular momentum lost to 
GWs and assuming e and r„ follow the Newtonian rela- 



tionship to energy and angular momentum. We approx- 
imate these losses as occurring in discrete steps during 
close-encounters. (This approximation will break down 
as e — >■ 0.) Trajectories computed with the NQA am- 
plitude and eccentricity-dependence [13, HI, [zH are also 
shown in the figure. The NQA execute many more orbits 
before merger. This clear departure from the NQA pre- 
diction at small rp due to strong field- GR effects should 
thus be apparent in the gravitational waveform. This 
model can also be used to predict approximately the 
critical-initial impact parameter for the BH-NS system 
to merge on the second close-encounter, the third close- 
encounter, and so on. We expect interesting dynamics 
around each of the thresholds up until the point where 
the system has undergone enough close-encounters that 
it circularizes. 



VI. CONCLUSION 

We have performed an initial survey of eccentric BH- 
NS mergers including the effects of black hole spin and 
varying the NS EOS. Though the limited number of val- 
ues we considered in this work does not begin to exhaust 
the parameter space, what is immediately apparent is the 
strong diversity in the resulting gravitational and matter 
dynamics. Though we have not yet studied the conse- 
quences of this on gravitational-wave detectability and 
parameter extraction, or possible electromagnetic coun- 
terparts, it is clear that the outcome can depend sensi- 
tively on the binary parameters and matter EOS. 

In order to understand the effects of large eccentricity, 
these simulations can be compared to the quasicircular 
BH-NS merger simulations of [H, which used the 
same piecewise polytrope equations of state. In [4^ . it 
was found that for a nonspinning BH and a 3:l-mass ra- 
tio, the resulting disk masses where 0.044, 0.0015, and 
< 10"^ Mq for the 2H, HB, and B EOSs, respectively. 
Hence eccentric mergers with certain ranges of impact 
parameters is one way to achieve significantly larger ac- 
cretion disks compared to the quasicircular case without 
BH spin. The effects of BH spin were considered in [ssj , 
where for a 4:l-mass ratio and the HB EOS they found 
disk masses of 0.024 and 0.18 Mq for BH spins of a = 0.5 
and 0.75, respectively (and even larger values for stifFer 
EOSs or lower mass ratios). They also only found a non- 
negligible amount of unbound material (> O.OIM0) only 
for stiff EOSs like the 2H. Thus, an important character- 
istic of eccentric mergers is the larger amount of ejected 
material found for some parameters. This could be a 
significant source of r-process elements (sol - fs^ , and give 
rise to EM counterparts, e.g. through nuclear decay of 
the radioactive r-process isotopes d, [sl] . The amount of 
energy radiated away in gravitational waves for the qua- 
sicircular case is somewhat higher than the values found 
here (for example, Eqw = 1.7% of the total mass for 
a 4:1 mass ratio, HB EOS, and a = 0.5 in [11]). The 
gravitational-wave signal is of course completely differ- 
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FIG. 13: The evolution of the eccentricity and periapsis separation of various 4 : 1 mass ratio BH-NS binaries that begin 
marginally unbound and undergo a series of close-encounters (large red points) before merging. For comparison, we also plot 
the results using the NQA expressions from [l^. [65j (magenta x's). We also plot the critical eccentricity for a given Vp for a 
close-encounter to result in merger (blue, dotted- line). Hence, the points below this curve correspond to merger events. From 
left to right the plots correspond to an initial BH-spin of a = —0.5, 0, and 0.5, respectively. 



ent for eccentric mergers, which are dominated by bursts 
from close-encounters or merger of the two compact ob- 
jects. 

Whether dynamical capture BH-NS mergers occur 
with sufficient frequency in the Universe to constitute a 
decent event rate for ground-based GW detectors is an- 
other question; at the least, failure to observe such events 
will place constraints on these sources, while if several are 
observed given the sensitivity of the outcome to proper- 
ties of the binary, they could be ideal environments to re- 
veal the structure of neutron stars. And again, we should 
emphasize that excessive fine tuning of parameters is not 
required for significant variability. Taking the relative 
velocity at large separations in a typical nuclear-cluster 
cusp as an example (w 1000 km/s), roughly 25% of such 
encounters will have < 10, (following the same line 
of reasoning discussed in the introduction for the per- 
centage of direct collisions) corresponding to the cases 
studied here. Certainly some of the most extreme exam- 
ples of mass transfer, large accretion disks, or multiple 
whirl orbits will be rare; perusing Tables UlTVl will give 
some idea of the distribution with rp. (Recall that the 
cross section scales linearly with Tp due to gravitational 
focusing, and that one would also expect some cases with 
larger initial rp than we followed through merger to ex- 
hibit similar variability.) It is not trivial to calculate a 
transition value of initial Tp above which the qualitative 
behavior at late times is described by a quasicircular in- 
spiral, though our study suggests at least a quarter of 
dynamical-capture binaries will merge with high eccen- 
tricity. 

For future work, we plan to investigate the detectabil- 
ity of the GW signals from BH-NS mergers that arise 
from dynamical capture in the strong-field regime. This 
will complement the first study of such systems pre- 
sented in [2l|, in that we intend to focus on the 



later stages of high-cccentricity mergers, including the 
merger/ringdown part of the signal. Future work also 
includes expanding the parameter space to different BH 
and NS masses, BH-spin orientations, and (as computa- 
tional resources permit) evolution of systems that exhibit 
more than two close-encounters before merger. Perform- 
ing higher-resolutions simulations will also be important 
to coming up with a more quantitatively accurate model 
of the behavior of the BH-NS binary near the threshold 
for merger during a given close-encounter. Doing so will 
not only require additional computational resources but 
a better method for creating initial data describing such 
eccentric binaries, since the superposition method used 
here imposes an effective floor on the accuracy that be 
reached. It would also certainly be interesting to investi- 
gate EM counterparts to these events; such simulations 
would require extensions to the code used here beyond 
the present GRHD. 



Acknowledgments 

We thank Adam Burrows, John Friedmann, Roman 
Gold, Janna Levin, Charalampos Markakis, Benjamin 
Lackey, Sean McWilliams and Richard O'Shaughnessy 
for useful conversations. This research was supported by 
the NSF through TcraGrid resources provided by NIGS 
under Grant No. TG-PHY100053, the Bradley Program 
fellowship (BGS), the NSF Graduate Research Fellow- 
ship under Grant No. DGE-0646086 (WE), NSF Grant 
No. PHY-0745779 (FP) and Grant No. PHY-1001515 
(BCS), and the Alfred P. Sloan Foundation (FP). Simu- 
lations were also run on the Woodhen cluster at Prince- 
ton University. 



17 



[1] A. Abramovici et al. Ligo: The laser interferometer grav- 
itational wave observatory. Science, 256:325, 1992. 

[2] Ramesli Narayan, Bohdan Paczynski, and Tsvi Piran. 
Gamma-ray bursts as the death throes of massive binary 
stars. Astrophys. J. , 395:L83, 1992. astro-ph/9204001. 

[3] N. Gehrels et al. A short 7-ray burst apparently asso- 
ciated with an elliptical galaxy at redshift z = 0.225. 
Nature (London) , 437:851-854, October 2005. 

[4] W. H. Lee, E. Ramirez-Ruiz, and J. Granot. A Compact 
Binary Merger Model for the Short, Hard GRB 050509b. 
Astrophys. J. Lett, 630:L165-L168, September 2005. 

[5] A. et al. Ran. Exploring the Optical Transient Sky with 
the Palomar Transient Factory. PASP, 121:1334-1351, 
December 2009. 

[6] N. Kaiser. Pan-STARRS: a wide-field optical survey 
telescope array. In J. M. Oschmann Jr., editor. Soci- 
ety of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, volume 5489 of Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference 
Series, pages 11-22, October 2004. 

[7] LSST Science Collaborations, P. A. AbeU, J. Alhson, 
S. F. Anderson, J. R. Andrew, J. R. P. Angel, L. Ar- 
mus, D. Arnett, S. J. Asztalos, T. S. Axelrod, and et al. 
LSST Science Book, Version 2.0. ArXiv e-prints, Decem- 
ber 2009. 

[8] B. D. Metzger and E. Berger. What is the Most Promis- 
ing Electromagnetic Counterpart of a Neutron Star Bi- 
nary Merger? ArXiv e-prints, August 2011. 
[9] N. Dalai, D. E. Holz, S. A. Hughes, and B. Jain. Short 
GRB and binary black hole standard sirens as a probe of 
dark energy. Phys. Rev. D , 74(6):063006, 2006. 

[10] Branson C. Stephens, William E. East, and Frans Preto- 
rius. Eccentric Black Hole-Neutron Star Mergers. Astro- 
phys. J. Lett., 737(1):L5, 2011. 

[11] R. Gold, S. Bernuzzi, M. Thierfelder, B. Bruegmann, 
and F. Pretorius. Eccentric binary neutron star merg- 
ers. ArXiv e-prints, September 2011. 

[12] Masaru Shibata and Keisuke Taniguchi. Coalescence of 
black hole-neutron star binaries. Living Reviews in Rel- 
ativity, 14(6), 2011. 

[13] W. H. Press and S. A. Teukolsky. On formation of 
close binaries by two-body tidal capture. Astrophys. J. , 
213:183-192, April 1977. 

[14] P.C. Peters and J. Mathews. Gravitational Radiation 
from Point Masses in a Keplerian Orbit. Phys. Rev., 
131:435-440, 1963. 

[15] R. M. O'Leary, B. Kocsis, and A. Loeb. Gravitational 
waves from scattering of stellar-mass black holes in galac- 
tic nuclei. Mon. Not. Roy. Astron. Soc, 395:2127-2146, 
June 2009. 

[16] W. H. Lee, E. Ramirez- Ruiz, and G. van de Ven. Short 
Gamma-ray Bursts from Dynamically Assembled Com- 
pact Binaries in Globular Clusters: Pathways, Rates, Hy- 
drodynamics, and Cosmological Setting. Astrophys. J. , 
720:953-975, September 2010. 

[17] J. Grindlay, S. Portegies Zwart, and S. McMillan. Short 
gamma-ray bursts from binary neutron star mergers in 
globular clusters. Nature Physics, 2:116-119, February 
2006. 

[18] R. Salvaterra, A. Cerutti, G. Chincarini, M. Colpi, 
C. Guidorzi, and P. Romano. Short Gamma-ray bursts: 



a bimodal origin? Mon. Not. Roy. Astron. Soc, 388:L6- 
L9, July 2008. 

[19] D. Guetta and L. Stella. Short 7-ray bursts and gravita- 
tional waves from dynamically formed merging binaries. 
Astron. and Astrophys., 498:329-333, May 2009. 

[20] K. Belczynski, M. Dominik, T. Bulik, R. O'Shaughnessy, 
C. Fryer, and D. E. Holz. The Effect of Metallicity on the 
Detection Prospects for Gravitational Waves. Astrophys. 
J. Lett, 715:L138-L141, June 2010. 

[21] B. Kocsis and J. Levin. Repeated Bursts from Relativis- 
tic Scattering of Compact Objects in Galactic Nuclei. 
ArXiv e-prints, September 2011. 

[22] J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy, 
T. Accadia, F. Acernese, C. Adams, R. Adhikari, 
P. Ajith, B. Allen, and et al. TOPICAL REVIEW: Pre- 
dictions for the rates of compact binary coalescences ob- 
servable by ground-based gravitational-wave detectors. 
Classical and Quantum Gravity, 27(17):173001, Septem- 
ber 2010. 

[23] E. Berger. A Short Gamma-ray Burst "No- host" Prob- 
lem? Investigating Large Progenitor Offsets for Short 
GRBs with Optical Afterglows. Astrophys. J. , 722:1946- 
1961, October 2010. 

[24] R. P. Church, A. J. Levan, M. B. Davies, and N. Tan- 
vir. Implications for the origin of short gamma-ray bursts 
from their observed positions around their host galaxies. 
Mon. Not. Roy. Astron. Soc, pages 355— 1-, March 2011. 

[25] J. D. M. Dewi, P. Podsiadlowski, and O. R. Pols. The 
spin period-eccentricity relation of double neutron stars: 
evidence for weak supernova kicks? Mon. Not. Roy. As- 
tron. Soc, 363:L71-L75, October 2005. 

[26] J. P. Norris, N. Gehrels, and J. D. Scargle. Heterogeneity 
in Short Gamma-Ray Bursts. Astrophys. J. , 735:23, July 
2011. 

[27] W. F. Domainko. Finding short GRB remnants in glob- 
ular clusters: the VHE gamma-ray source in Terzan 5. 
Astron. and Astrophys., 533:L5, September 2011. 

[28] T. A. Thompson. Accelerating Compact Object Merg- 
ers in Triple Systems with the Kozai Resonance: A 
Mechanism for "Prompt" Type la Supernovae, Gamma- 
Ray Bursts, and Other Exotica. Astrophys. J. , 741:82, 
November 2011. 

[29] Linqing Wen. On the Eccentricity Distribution of Coa- 
lescing Black Hole Binaries Driven by the Kozai Mech- 
anism in Globular Clusters. Astrophys. J., 598:419-430, 
2003. 

[30] J. M. Lattimer and D. N. Schramm. Black-hole-neutron- 
star collisions. Astrophys. J. Lett., 192:L145-L147, 
September 1974. 

[31] S. Rosswog, F.K. Thielemann, M.B. Davies, W. Benz, 
and T. Piran. Coalescing Neutron Stars: A Solution to 
the R-Process Problem? 1998. 

[32] Li-Xin Li and Bohdan Paczynski. Transient events from 
neutron star mergers. Astrophys. J., 507:L59, 1998. 

[33] B. D. Metzger, G. Martmez-Pinedo, S. Darbha, 
E. Quataert, A. Arcones, D. Kasen, R. Thomas, P. Nu- 
gent, I. V. Panov, and N. T. Zinner. Electromagnetic 
counterparts of compact object mergers powered by the 
radioactive decay of r-process nuclei. Mon. Not. Roy. 
Astron. Soc, 406:2650-2662, August 2010. 

[34] Charles F. Gammie, Stuart L. Shapiro, and Jonathan C. 



18 



McKinney. Black Hole Spin Evolution. Astrophys. J., [51 
602:312-319, 2004. 

[35] J. M. Miller, M. C. Miller, and C. S. Reynolds. The 

Angular Momenta of Neutron Stars and Black Holes as [52 
a Window on Supernovae. ArXiv e-prints, 2011. 

[36] Z. B. Etienne, Y. T. Liu, S. L. Shapiro, and T. W. [53 
Baumgarte. General relativistic simulations of black- 
hole-neutron-star mergers: Effects of black-hole spin. 
Phys. Rev. D , 79(4):044024, 2009. [54 

[37] Francois Foucart, Matthew D. Duez, Lawrence E. Kid- 
der, and Saul A. Teukolsky. Black hole-neutron star 
mergers: effects of the orientation of the black hole spin. 
Phys. Rev., D83:024005, 2011. [55 

[38] Koutarou Kyutoku, Hirotada Okawa, Masaru Shibata, 

and Keisuke Taniguchi. Gravitational waves from spin- [56 
ning black hole-neutron star binaries: dependence on 
black hole spins and on neutron star equations of state. 
Phys. Rev., D84:064018, 2011. [57 

[39] Sarvnipun Chawla et al. Mergers of Magnetized Neutron 
Stars with Spinning Black Holes: Disruption, Accretion 
and Fallback. Phys. Rev. Lett, 105:111101, 2010. [58 

[40] J. S. Read, C. Markakis, M. Shibata, K. Uryii, J. D. E. 

Creighton, and J. L. Friedman. Measuring the neutron [59 
star equation of state with gravitational wave observa- 
tions. Phys. Rev. D , 79(12):124033, 2009. 

[41] M. Shibata and K. Uryti. Merger of black hole-neutron 

star binaries: Nonspinning black hole case. Phys. Rev. D [60 
, 74(12):121503, 2006. 

[42] Matthew D. Duez, Francois Foucart, Lawrence E. Kid- 
der, Christian D. Ott, and Saul A. Teukolsky. Equation [61 
of state effects in black hole-neutron star mergers. Class. 
Quant. Grav., 27:114106, 2010. 

[43] Francesco Pannarale, Luciano RezzoUa, Frank Ohme, [62 
and Jocelyn S. Read. Will black hole-neutron star bi- 
nary inspirals tell us about the neutron star equation of 
state? ArXiv e-prints, 2011. 

[44] K. Kyutoku, M. Shibata, and K. Taniguchi. Gravita- [63 
tional waves from nonspinning black hole-neutron star 
binaries: Dependence on equations of state. Phys. Rev. [64 
D , 82(4):044049— August 2010. 

[45] Benjamin D. Lackey, Koutarou Kyutoku, Masaru Shi- 
bata, Patrick R. Brady, and John L. Friedman. Extract- [65 
ing equation of state parameters from black hole- neutron 
star mergers. I. Nonspinning black holes. ArXiv e-prints, 
2011. [66 

[46] S. Chawla, M. Anderson, M. Besselman, L. Lehner, S. L. 
Liebling, P. M. Motl, and D. Neilsen. Mergers of Magne- 
tized Neutron Stars with Spinning Black Holes: Disrup- 
tion, Accretion and Fallback. ArXiv e-prints, 2010. [67 

[47] Z. B. Etienne, Y. T. Liu, V. Paschalidis, and S. L. 
Shapiro. General relativistic simulations of black hole- 
neutron star mergers: Effects of magnetic fields. ArXiv 
e-prints, December 2011. [68 

[48] F. Foucart, M. D. Duez, L. E. Kidder, M. A. Scheel, 
B. Szilagyi, and S. A. Teukolsky. Black hole-neutron star 
mergers for 10 solar mass black holes. ArXiv e-prints, [69 
November 2011. 

[49] W. E. East, F. Pretorius, and B. C. Stephens. Hydrody- 
namics in full general relativity with conservative AMR. [70 
ArXiv e-prints, December 2011. 

[50] F. Lofiler, L. RezzoUa, and M. Ansorg. Numerical 

evolutions of a black hole-neutron star system in full [71 
general relativity: Head-on collision. Phys. Rev. D , 
74(10):104018, 2006. 



H. Friedrich. On the hyperbolicity of Einstein's and other 
gauge field equations. Communications in Mathematical 
Physics, 100:525-543, December 1985. 
D. Garfinkle. Harmonic coordinate method for simulating 
generic singularities. Phys. Rev. D , 65(4):044029, 2002. 
F. Pretorius. Numerical relativity using a general- 
ized harmonic decomposition. Class, and Quant. Grav., 
22:425, 2005. 

Carsten Gundlach, Jose M. Martin-Garcia, Gioel Cal- 
abrese, and Ian Hinder. Constraint damping in the Z4 
formulation and harmonic gauge. Class. Quant. Grav., 
22:3767-3774, 2005. 

F. Pretorius. Evolution of Binary Black-Hole Spacetimes. 
Phys. Rev. Lett. , 95(12):121101, 2005. 

Matthew W. Choptuik and Frans Pretorius. Ultra Rela- 
tivistic Particle Collisions. Phys. Rev. Lett., 104:111101, 
2010. 

Lee Lindblom and Bela Szilagyi. An Improved Gauge 
Driver for the GH Einstein System. Phys. Rev., 
D80:084019, 2009. 

A. Harten, P.D. Lax, and B.J. van Leer. SI AM Rev., 
25:35, 1983. 

A. Tchekhovskoy, J. C. McKinney, and R. Narayan. 
WHAM: a WENO-based general relativistic numerical 
scheme - 1. Hydrodynamics. Man. Not. Roy. Astron. Soc., 
379:469, 2007. 

Marsha J Berger and Joseph Oliger. Adaptive mesh 
refinement for hyperbolic partial differential equations. 
Journal of Computational Physics, 53(3) :484, 1984. 
M.J. Berger and P. Colella. Local adaptive mesh refine- 
ment for shock hydrodynamics. .Journal of Computa- 
tional Physics, 82(1):64, 1989. 

G. D. Quinlan and S. L. Shapiro. The collapse of dense 
star clusters to supermassive black holes - Binaries and 
gravitational radiation. Astrophys. J. , 321:199-210, Oc- 
tober 1987. 

James M. Lattimer and M. Prakash. What a Two Solar 
Mass Neutron Star Really Means. ArXiv e-prints, 2010. 
Feryal Ozel, Dimitrios Psaltis, Ramesh Narayan, and Jef- 
frey E. McChntock. The Black Hole Mass Distribution 
in the Galaxy. ArXiv e-pnnts, 2010. 
M. Turner. Gravitational radiation from point-masses 
in unbound orbits - Newtonian results. Astrophys. J. , 
216:610-619, 1977. 

Paul Demorest, Tim Pennucci, Scott Ransom, Mallory 
Roberts, and Jason Hessels. Shapiro delay measurement 
of a two solar mass neutron star. Nature, 467:1081-1083, 
2010. 

Kostas Glampedakis and Daniel Kennefick. Zoom and 
whirl: Eccentric equatorial orbits around spinning black 
holes and their evolution under gravitational radiation 
reaction. Phys. Rev., D66:044002, 2002. 
Frans Pretorius and Deepak Khurana. Black hole merg- 
ers and unstable circular orbits. Class. Quant. Grav., 
24:S83-S108, 2007. 

M. J. Rees. Tidal disruption of stars by black holes of 10 
to the 6th-10 to the 8th solar masses in nearby galaxies. 
Nature (London) , 333:523-528, June 1988. 
W. Kastaun, B. Willburger, and K. D. Kokkotas. Satu- 
ration amplitude of the f-mode instability. Phys. Rev. D 
, 82(10):104036, November 2010. 

K.D. Kokkotas, T.A. Apostolatos, and N. Andersson. 
The Inverse problem for pulsating neutron stars: A 'Fin- 
gerprint analysis' for the supranuclear equation of state. 



19 



Mon.Not.Roy.Astron.Soc, 320:307-315, 2001. 

[72] Christopher P. L. Berry and Jonathan R. Gair. Gravi- 
tational wave energy spectrum of a parabolic encounter. 
Phys. Rev., D82:107501, 2010. 

[73] K. Glampedakis and D. Kennefick. Zoom and whirl: 
Eccentric equatorial orbits around spinning black holes 



and their evolution under gravitational radiation reac- 
tion. Phys. Rev. D , 66(4):044002, 2002. 
[74] Neil J. Cornish and Janna J. Levin. Lyapunov timescales 
and black hole binaries. Class. Quant. Grav., 20:1649- 
1660, 2003. 



